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Abstract 

We construct new black hole solutions in Einstein- Yang-Mills theory. They are 
static, axially symmetric and asymptotically flat. They are characterized by their 
horizon radius and a pair of integers (k,n), where k is related to the polar angle 
and n to the azimuthal angle. The known spherically and axially symmetric EYM 
black holes have k = 1. For k > 1, pairs of new black hole solutions appear above a 
minimal value of n, that increases with k. Emerging from globally regular solutions, 
they form two branches, which merge and end at a maximal value of the horizon 
radius. The difference of their mass and their horizon mass equals the mass of the 
corresponding regular solution, as expected from the isolated horizon framework. 
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1 Introduction 



The well-known regular Bartnik-McKinnon (BM) solutions pQ and the corresponding non- 
Abelian black hole solutions [2j, are asymptotically flat, static spherically symmetric solu- 
tions of SU(2) Einstein- Yang-Mills (EYM) theory. They are unstable solutions, sphalerons 

and are characterized by the number of nodes of the gauge field. Besides these spheri- 
cally symmetric solutions there are also asymptotically flat, static regular and black hole 
solutions, which possess only axial symmetry [3]. These are characterized by two integers, 
the node number of their gauge field function(s), and the winding number with respect to 
the azimuthal angle, denoted n. The spherically symmetric solutions have winding number 
n — 1, while the only axially symmetric solutions have winding number n > 1. The n > 1 
black hole solutions possess an event horizon with a slight elongation along the symmetry 
axis Hj. All these EYM black hole solut ions exist for arbitrarily large horizon size. 

In SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory with a triplet Higgs field, gravi- 
tating monopole solutions and black holes with monopole hair arise El • As in EYM 
theory, the regular EYMH solutions are characterized by two integers, the node number of 
the gauge field and the azimuthal winding number n, which corresponds to the topological 
charge of the monopoles. The black hole solutions are characterized in addition by their 
horizon size. These EYMH black hole solutions exist only up to a maximal value of the 
horizon size [El IE], unlike the known EYM black hole solutions. 

EYMH theory allows for further static axially symmetric solutions, representing grav- 
itating monopole-antimonopole pair, chain and vortex solutions pH El E], which can be 
characterized by the azimuthal winding number n, and by a second integer m, related to 
the polar angle. For the monopole-antimonopole chains, which arise in flat space for n = 1 
and 2, the integer m corresponds to the number of nodes of the Higgs field (and thus the 
number of poles on the symmetry axis), while in vortex solutions, which arise in flat space 
for winding number n > 2, the Higgs field vanishes (for even m) on m/2 rings centered 
around the symmetry axis jHj- To all these regular solutions associated black hole solutions 
should exist, obtained so far only for a monopole-antimonopole pair [TO] . 

In the limit of vanishing Higgs expectation value, EYMH solutions approach (after 
rescaling) EYM solutions |HJ . Interestingly, when n > 4 and m > 4, new regular EYM 
solutions appear as limiting solutions E]. These static axially symmetric solutions 
have been characterized by the integers (k, n), where 2k = m. They have been constructed 
numerically for k = 2, n > 4, and k = 3, n > 6 [llj. Unlike the k = 1 EYM solutions, 
the k > 1 solutions always appear in pairs. In this letter we construct the corresponding 
EYM black hole solutions. As expected, a branch of black hole solutions is associated with 
each regular solution. Intriguingly, the two branches of black hole solutions from a pair of 
regular solutions merge and end at a maximal value of the horizon size. 

In the isolated horizon framework ^2j non-Abelian black hole solutions can 

be interpreted as bound states of regular solutions and Schwarzschild black holes [Tlj . 
Furthermore, the isolated horizon framework yields a relation for the mass of non-Abelian 
black hole solutions, representing it as the sum of the mass of the regular solution and the 
horizon mass of the black hole solutions [T^]. We here show that this relation is also valid 
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for the new black holes. In particular, the two regular solutions of a given k and n are 
connected via this mass formula. 

In section II we present the EYM action, the axially symmetric ansatz and the boundary 
conditions. In section III we address their asymptotic and horizon properties. We present 
our numericl results in section IV, and we give our conclusions in section V. 



2 Action and Ansatz 

We consider the SU(2) EYM action 

with Ricci scalar R, field strength tensor 

F^ = d ll A v -d v A li + ie[A lt ,A v \ , (2) 

gauge potential A^ = r a A^/2, and gravitational and Yang-Mills coupling constants G and 
e, respectively. Variation of the action with respect to the metric g^ v leads to the 
Einstein equations, variation with respect to the gauge potential A^ to the gauge field 
equations. 

In isotropic coordinates the static axially symmetric metric reads 0] 

ds > = - fd e + ^ + ^fde* + , (3) 

where the metric functions /, m and I are functions of the coordinates r and 9, only. The 
z-axis (9 = 0, 7r) represents the symmetry axis. Regularity on the z-axis requires m = I 
there. 

For the gauge field we employ the ansatz [U IHJ EO HI] 



AJx^ = — 
2er 



r£ (H^r + (1 - H 2 ) rd9) - n (r?> k H 3 + T^ k H 4 ) r sin 9d<f)] . (4) 



Here the symbols r" ,fc , r^' fc and denote the dot products of the cartesian vector of Pauli 



matrices, r = (r x , r y , r z ), with the spatial unit vectors 

e r n,k = (sin k9 cosmp, sin k9 slump, cos k6) , 
e 9 n ' h = (cos kQ cos rup, cos k9 smrup, — sin k9) , 
ej 1 = (— sinn<£>, cosmp, 0) , (5) 

respectively. The gauge field functions ifj, % — 1 — 4, depend on the coordinates r and 9, 
only. For k = n = 1 and H\ = H% = 0, H2 = 1 — H4 = w(r) the spherically symmetric 
black hole solutions |2j are recovered, while for k = 1, n > 1, one obtains the axially 
symmetric solutions of 0]. The new black hole solutions reported here are obtained for 
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k > 1. In the limit of vanishing horizon radius they converge pointwise to the globally 
regular solutions with the same integers k and n jTT] . The globally regular solutions are 
related to EYMH solutions with m = 2k in the limit of vanishing Higgs field [TT| 19]. 
The ansatz is form-invariant under the Abelian gauge transformation [1] 

U = exp Q^r(r, 0)) . (6) 

We fix the gauge by choosing the gauge condition jU |H1 El HI] 

rd r H t - d e H 2 = . (7) 

To obtain asymptotically flat solutions which are regular at the horizon and possess 
the proper symmetries, we need to impose appropriate boundary conditions jU El 13 HI] • 
The horizon of the non- Abelian black hole solutions resides at a surface of constant radial 
coordinate r = th- At the horizon we impose the boundary conditions 

/ = m = I = o , Fi = 0, d r H 2 = d r H 3 = d r H A = , (8) 

at infinity we impose 

/ = m = l = 1 , H x = #3 = , H 2 = l-2k , H 4 = 2 sm(k9) / sin 9 , (9) 

and on the z-axis we impose 

d e f = d e m = d e l = 0, H 1= H 3 = 0, d e H 2 = d e H A = . (10) 



3 Properties 



We introduce the dimensionless coordinate x = ,? „ r and horizon radius x-h = ,f ^ r-u. 
Defining the mass M of the black hole solutions via the Komar integral, the dimensionless 
mass n = ^|=jM, is determined by the derivative of the metric function / at infinity 

^=^i™x 2 ^/. (11) 

From the equations of motion it follows [3] , that the Kretschmann scalar is finite at the 
horizon, and that the surface gravity k [15 , 

k 2 = -(l/4)sV'(^ tt )(%tt) , (12) 

is constant, as required by the zeroth law of black hole physics. Expansion of the metric 
functions near the horizon in the form 

f(x, 9) = x 2 f 2 (9) + O (5 3 ) , m(x, 9) = £ 2 m 2 (9) + O [x 3 
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where x — (x/x-r — 1), yields the dimensionless surface gravity k 



h{0) 



related to 



x^Jm 2 {9) 



k by k = ke/yAftG. 

We introduce the area parameter x& ITS] , defined via the dimensionless area of the 
black hole horizon A, 



A = 2ttx 



ii 



d6 sin 9 



f 



47TX1 



(13) 



The deformation of the horizon is revealed, when the circumference of the horizon along 
the equator, L e , is compared to the circumference of the horizon along the poles, L p , 



2- 



dip 



-x sin 9 



Lp — 2 



x=xn,8=n/2 



de 



-x 



(14) 



x=xn,tp=const. 



since the black hole solutions have L p ^ L e (in general). 

The isolated horizon framework ^21 EH] yields an intriguing relation between the ADM 
mass /i of a black hole with area parameter xa and the mass /i reg of the corresponding 
globally regular solution [P2] . 

yU = flA + yWreg , (15) 

where the (dimensionless) horizon mass //a is defined via 



th^x^jx^dx^ . 



(16) 



The isolated horizon formalism further suggests to interpret a non-Abelian black hole 
as a bound state of a regular solution and a Schwarzschild black hole 



/i — /i reg + /is + /ibind 



(17) 



where /is = ^a/2 is the ADM mass of the Schwarzschild black hole with area parameter 
xa, and /ibmd represents the binding energy of the system, 



/ibind — A 4 a — A*s • 



4 Results 

Subject to the above boundary conditions, we solve the system of seven coupled non-linear 
partial differential equations numerically. To map spatial infinity to the finite value x = 1, 
we employ the radial coordinate 

- x-x u 

x = — . (19) 

1 + x y ' 

The numerical calculations are based on the Newton-Raphson method, and are performed 
with help of the program FIDISOL ^H]- The equations are discretized on a non-equidistant 
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grid in x and 6. Typical grids used have sizes 70 x 30, covering the integration region 
0<x<land0<#< tt/2. 

We construct numerically the black hole solutions for (A; = 1, n — 1 — 8), (k — 2, 
n = 4 — 8) and (k — 3, n = 6 — 8). For k — 1 and fixed n a branch of black hole solutions 
emerges from the single globally regular solution, when the horizon radius ih is increased 
from zero jl] . This branch of black hole solutions exists for arbitrarily large horizon size. 
When k > 2 and n > 2k, however, a pair of globally regular solutions exists for a given set 
(k, n). Thus we find two branches of black hole solutions emerging from the corresponding 
globally regular solutions, when the horizon radius is increased from zero. Surprisingly, the 
two branches do not exist for arbitrarily large horizon size. Instead they merge and end 
at a maximal value of the horizon size, i.e., at the maximal value of the area parameter 
%A,max, which depends on k and n. 

We show the ADM mass /i as a function of the area parameter xa in Fig. 1. The 
mass increases monotonically with increasing xa- For the new solutions with k > 2, the 
maximal horizon size increases with increasing n. At the same time, the mass difference 
of the two regular solutions increases, indicating a possible correlation between the soliton 
mass difference and the maximal horizon size. 

Considering the shape of the horizon, we observe that the deviation from spherical 
symmetry is small, though. We observe a ratio of circumferences for k = 1 solutions of 
up to L e /L p = 0.988, and for k = 2 and k = 3 solutions of up to L e /L p = 0.984 and 
L e /L p = 0.979, respectively, as seen in Fig. 2. Thus the new black holes have a slightly 
more deformed horizon. 

The inverse surface gravity 1/k of the black hole solutions is exhibited in Fig. 3. For 
k = 1 black hole solutions the inverse surface gravity increases monotonically (except for 
n = 1) as a function of the area parameter xa- For k = 2 and k = 3 black holes it increases 
along the lower branch, and beyond the transition to the upper branch it decreases back 
to zero. 

Our numerical results indicate, that the mass relation (fT3j) . obtained in the isolated 
horizon formalism, also holds for the new non-Abelian black hole solutions. (In Fig. 1 for 
each branch the corresponding regular solution is the reference point for the integration.) 
Consequently, the mass of the regular solution on the upper branch is related to the mass 
of the regular solution on the lower branch via the horizon mass integral, performed along 
both branches. A similar result was obtained previously for EYMH black holes with dipole 
hair JTUI - 

In [Hj the mass formula /z reg = 1/2 J^°(l — xak)cIxa was derived and shown to hold 
for spherically symmetric EYM solutions in . We checked that the same formula also 
holds for the k = 1 axially symmetric solutions. 

In Fig. 4 we illustrate the binding energy of the non-Abelian black hole solutions. The 
binding energy is always calculated w.r.t. the corresponding regular solution. Consequently, 
for k = 2 and k = 3 black holes the binding energy along the upper branches is smaller 
than the binding energy along the lower branches. The difference in binding energy of the 
black hole solutions for a given set of (k, n) at the maximal value of the area parameter 
represents the mass difference of the corresponding pair of regular solutions. 



6 






Fig. 1 (a) The ADM mass \x is shown as function of the area parameter xa for k — 1, 
n = 1— 8 (symbols x). Also shown is /XA+/x reg (lines), (b) Same as (a) for k = 2, n = 4 — 8. 
(c) Same as (a) for k = 3, n = 6 — 8. Fig. 2 Same as Fig. 1 for the ratio L e /L p . 
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5 Conclusions 

We have constructed numerically new static axially symmetric black hole solutions of EYM 
theory and investigated their properties. These solutions are characterized by the set of 
integers (k,n), related to the polar and azimuthal angles, respectively. In particular, we 
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have obtained black hole solutions for k = 2, n = 4 — 8, and k = 3, n = 6 — 8. 

The (l,n) solutions (with one node) exist most likely for any integer n > 1 and for 
arbitrary horizon size. The new (2, n) and (3, n) solutions have lower bounds on n, n = 2k, 
and upper bounds on the horizon size. For each allowed set (k, n) two branches of black 
hole solutions emerge from the two globally regular solutions and merge and end at the 
maximal value of the horizon size. The existence of an upper bound of the horizon radius 
is a surprising and interesting new feature for EYM solutions. Previously a maximal value 
of the horizon size was known only in theories with an in-built length scale already in 
Minkowski space, hence it was conjectured, that "there is no bound on the horizon radius 
of hairy, static EYM black holes" [T2*] . The underlying reason for the occurrence of this 
maximal horizon size for the new EYM black holes is yet to be understood. 

We expect, that the (2, n) and (3, n) regular and black hole solutions represent only the 
first sequences of new solutions, and conjecture the existence of (k, n) regular and black 
hole solutions also for higher values of k. The (sperically symmetric) k — 1 black holes are 
unstable and there is all reason to believe, that the new (k,n) solutions are unstable 
as well. 

Considering the new static axially symmetric solutions from the isolated horizon for- 
malism point of view, we have verified the mass relation between the regular and the black 
hole solutions, showing that the black hole mass is given by the sum of the mass of the 
regular solution and the horizon mass. In particular, the mass of the regular solution on 
the upper branch is related to the mass of the regular solution on the lower branch via 
the horizon mass integral, performed along both branches. Interpreting the non-Abelian 
black holes as bound states of regular solutions and Schwarzschild black holes [H], we have 
obtained the binding energy of these bound systems. 

We note that the globally regular solutions were first obtained in EYMH theory as 
limiting solutions for vanishing Higgs field [01 ITT] . Hence we expect that EYMH black 
holes can be found from the new EYM black holes by gradually switching on the Higgs 
field. 

Rotating EYM black hole solutions based on the static k = 1 black hole solutions are 
known ^H]- It appears straightforward to construct rotating EYM black hole solutions 
based on the static k = 2 and k = 3 black hole solutions. In contrast, rotating regular 
EYM solutions do not appear to exist JO] > though recently the first rotating regular EYMH 
solutions have been constructed |T[|1 • 
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